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Abstract 

We search for CP violating asymmetries in the charmless hadronic B 
meson decays to K^if?, K^tt , K^, K±rf , and wvr*. With the full CLEO 
II and CLEO II. V datasets statistical precision on „4cp is in the range of 
±0.12 to ±0.25 depending on the mode. All quoted results are preliminary. 
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1 Introduction 



CP violating phenomena arise in the Standard Model because of the single complex pa- 
rameter in the quark mixing matrix. Q Such phenomena are expected to occur widely 
in B meson decays and are the incentive for most of the current -B-physics initiatives 
in the world. As of yet there is little direct experimental evidence. CDF's recent 
determination^ of sin 2/3 at the 2a level, which followed upon earlier less sensitive 
searches by both CDF||] and OPALQ, is consistent with expectations for mixing- induced 
Standard Model CP violation and to date the only evidence of CP effects in B mesons. 

Direct CP violation however is also anticipated to play a prominent role in the CP 
phenomena of B decay. To date the only published search for direct CP violation in B 
decay is the recent CLEO limit || on Aqv in b — ► 57. CP asymmetries may show up 
in any mode where there are two or more participating diagrams which differ in weak 
and strong phases. B — * Kir modes, for instance, involve b —* u tree diagrams carrying 
the weak phase Aig(V* b V us ) = 7 and b — » s penguin diagrams carrying the weak phase 
Arg(T^Vfcs) = 7r. Though the branching ratios are small, such cases are experimentally 
straightforward to search for. Rate differences between B — ► K + ir~ and B — > K~ir + 
decays would be unambiguous signals of direct CP violation if seen. 

We report here five searches for direct CP violation in charmless hadronic B decay 
modes, based on the full CLEO II and CLEO II. V datasets which together comprise 
9.66 million BB events. The modes searched for are the three Kir modes, K^tt^, 
K^tt , KgTr^, the mode K rf, and the vector-pseudoscalar mode ujiv . In all but the 
first case the flavor of the parent b or b quark is tagged simply by the sign of the high 
momentum charged hadron; for K^tt^ one must further identify K and it. 

In what follows we will refer when needed to the generic final state from b and b as 
/ and / respectively. The corresponding event yields of signal and background we will 
label as S, S, B, and B. We define the sign of -4cp with the following convention: 

_ Brjp -> /) - Br(b -> f) S -S 
^ CP Br(b -> /) + Br(b - /) S + S {> 

The statistical precision one can achieve in a measurement of .4cp depends on the 
signal yield S = S + S, the CP-symmetric background B = B + B, any correlation \ 
that might exist between the measurements of / and /, and of course on Acv itself: 



^cp S \ 1 ^ S I 1 - ^ CP 



2 1-Acp 2 L , B (1 + Acp 




In most cases there is no chance of confusing / and / so x = 0. However for / = 
K~ir + and / = K + n~, a small degree of crossover is possible due to imperfect particle 
identification; we find % = —0.11 for this case. For % ~ 0) B/S ~ 1, and Ac? ~ one 
has an easy rule of thumb, a ~ y/2/S. For S ~ 100 this means one expects statistical 
precision in the neighborhood of ±0.15. As will be discussed later, systematic errors 
are small, and consequently the precision in .Acp measurements can be expected to be 
dominated by statistical errors for a long time to come. As can be seen in Eq. the 
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statistical error is in turn dominated by the leading l/y/S coefficient which reminds us 
that the only path to better Aqp measurements will be more data. Improvements to 
analysis technique that reduce B/S or x have less impact. 

2 Theoretical Expectations 

The existence of a CP violating rate asymmetry depends on having both two different 
CP nonconserving weak phases and two different CP conserving strong phases. The 
former may arise from either the Standard Model CKM matrix or from new physics, 
while the latter may arise from the absorptive part of a penguin diagram or from final 
state interaction effects. The difficulty of calculating strong interaction phases, particu- 
larly when long distance non-perturbative effects are involved, largely precludes reliable 
predictions. Under well-defined model assumptions, however, numerical estimates may 
be made and the dependence on both model parameters and CKM parameters can be 
probed. A recent and comprehensive review of CP asymmetries under the assumption of 
generalized factorization has been published by Ali et al Q. We quote in Table |] their 
predictions for the modes examined in this paper. 



Mode A C p 



B - 


■* 


0.037 ^0.106 


B - 




0.026 -> 0.092 


B - 


-> K ^ 


0.015 


B - 


-> K^rf 


0.020 -> 0.061 


B - 


-0.120 -> 0.024 



Table 1: CP asymmetry predictions. The range of ^4cp values reflects a range of 
model parameters. CKM parameters p and rj are set to p = 0.12, rj = 0.34. Taken 
from Ali et al, Ref. 6, with signs changed to match the convention of Eq. 1. 



If final state interactions are not neglected, however, strong phases as large as 90° are 
not ruled out and |-4cp| could reach 0.44 in favorable cases. [^] In this nonperturbative 
regime numerical predictivity is limited, but a variety of relationships among asymmetries 
or /, / rate differences can be found in the literature.^] 

3 Data Set, Detector, Event Selection 

The data set used in this analysis was collected with the CLEO II and CLEO II. V 
detectors at the Cornell Electron Storage Ring (CESR). It consists of 9.1 fb^ 1 taken at 
the T(4S) (on-resonance) and 4.5 fb _1 taken below BB threshold. The below-threshold 
sample is used for continuum background studies. The on-resonance sample contains 
9.66 million BB pairs. This is a factor 2.9 increase in the number of BB pairs over 



the published measurements of the modes considered here [15], [16 ],[!?]. In addition, 
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the CLEO II. V data set, which has significantly improved particle identification and 
momentum resolution as compared with CLEO II, now dominates the data set. 

CLEO II and CLEO II. V are general purpose solenoidal magnet detectors, described 
in detail elsewhere [pj. In CLEO II, the momenta of charged particles are measured in 
a tracking system consisting of a 6-layer straw tube chamber, a 10-layer precision drift 
chamber, and a 51-layer main drift chamber, all operating inside a 1.5 T superconducting 
solenoid. The main drift chamber also provides a measurement of the specific ionization 
loss, dE/dx, used for particle identification. For CLEO II. V the 6-layer straw tube 
chamber was replaced by a 3-layer, double-sided silicon vertex detector, and the gas in 
the main drift chamber was changed from an argon-ethane to a helium-propane mixture. 
Photons are detected using 7800-crystal CsI(Tl) electromagnetic calorimeter. Muons are 
identified using proportional counters placed at various depths in the steel return yoke 
of the magnet. 

Charged tracks are required to pass track quality cuts based on the average hit 
residual and the impact parameters in both the r — <fi and r — z planes. Candidate K s 
are selected from pairs of tracks forming well-measured displaced vertices. Furthermore, 
we require the Kg momentum vector to point back to the beam spot and the 7r + 7r~ 
invariant mass to be within 10 MeV, two standard deviations (<r), of the K s mass. 
Isolated showers with energies greater than 40 MeV in the central region of the Csl 
calorimeter and greater than 50 MeV elsewhere, are defined to be photons. Pairs of 
photons with an invariant mass within 2.5a of the nominal 7r° (77) mass are kinematically 
fitted with the mass constrained to the nominal ir° (77) mass. To reduce combinatoric 
backgrounds we require the lateral shapes of the showers to be consistent with those 
from photons. To suppress further low energy showers from charged particle interactions 
in the calorimeter we apply a shower-energy-dependent isolation cut. 

Charged particles are identified as kaons or pions using dE/dx. Electrons are rejected 
based on dE/dx and the ratio of the track momentum to the associated shower energy 
in the Csl calorimeter. We reject muons by requiring that the tracks do not penetrate 
the steel absorber to a depth greater than seven nuclear interaction lengths. We have 
studied the dE/dx separation between kaons and pions for momenta p ~ 2.6 GeV/c in 
data using D° — ► K~tt + (tt°) decays; we find a separation of (1.7 ± 0.1) a for CLEO II 
and (2.0 ± 0.1) a for CLEO II.V. 

Resonances are reconstructed through the decay channels: rj' — * t]tt + tt~ with 77 — * 77; 
7/ — > p7 with p — > 7T + 7r - ; and u> — > 7r + 7r~7r°. 

4 Analysis 

The *4cp analyses presented are intimately related to the corresponding branching ratio 
determinations presented in separate contributions to this conference. [|| We summarize 
here the main points of the analysis. 

We select hadronic events and impose efficient quality cuts on tracks, photons, ir° can- 
didates, and Kg candidates. We calculate a beam-constrained B mass M = yj E^ — p 2 B , 
where pb is the B candidate momentum and E^ is the beam energy. The resolution 
in M ranges from 2.5 to 3.0 MeV/c 2 , where the larger resolution corresponds to the 
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B ± -» /i ± vr° decay. We define AE = E x + E 2 - E h , where £1 and E 2 are the energies of 
the daughters of the B meson candidate. The resolution on AE is mode-dependent. For 
final states without photons the AE resolution for CLEO II.V(II) is 20(26) MeV. Most 
other modes are only slightly worse but for the — > h^ir analysis, the AE resolution 
is worse by about a factor of two and becomes asymmetric because of energy loss out 
of the back of the Csl crystals. The energy constraint also helps to distinguish between 
modes of the same topology. For example, AE for B —* K + tt~, calculated assuming 
B — > 7T + 7r _ , has a distribution that is centered at —42 MeV, giving a separation of 
2.1(1.6)<7 between B -> K+tt' and B -> vt+tt- for CLEO II.V(II). We accept events 
with M within 5.2 - 5.3 GeV/c 2 and \ AE\ < 200 MeV. The AE requirement is loosened 
to 300 MeV for the B ± -> h ± ir° analysis. This fiducial region includes the signal region, 
and a sideband for background determination. Similar regions are included around each 
of the resonance masses (r/, r/, and to) in the likelihood fit. For the 7/ — > /r/ case, the p 
mass is not included in the fit; we require 0.5 GeV < < 0.9 GeV. 

We have studied backgrounds from b — > c decays and other b —* u and b —* s decays 
and find that all are negligible for the analyses presented here. The main background 
arises from e + e~ — ► qq (where q = u,d,s,c). Such events typically exhibit a two-jet 
structure and can produce high momentum back-to-back tracks in the fiducial region. 
To reduce contamination from these events, we calculate the angle 8 sp h between the 
sphericity axis[ll| of the candidate tracks and showers and the sphericity axis of the rest 
of the event. The distribution of cos # sp h is strongly peaked at ±1 for qq events and 
is nearly flat for BB events. We require |cos# sp h| < 0.8 which eliminates 83% of the 
background. For 7/ and oj modes the cut is made at 0.9. 

Additional discrimination between signal and qq background is provided by a Fisher 
discriminant technique as described in detail in Ref. [13]. The Fisher discriminant is a 
linear combination T = J2iLi a iUi where the coefficients a, are chosen to maximize the 
separation between the signal and background Monte-Carlo samples. The 11 inputs, yi, 
are | cos6 ca nd\ (the cosine of the angle between the candidate sphericity axis and beam 
axis), the ratio of Fox- Wolfram moments H 2 /Hq [14|, and nine variables that measure 
the scalar sum of the momenta of tracks and showers from the rest of the event in nine 
angular bins, each of 10°, centered about the candidate's sphericity axis. For the rf and 
ui modes, | cos#b| (the angle between the B meson momentum and beam axis), is used 
instead of H 2 /-Ho- 
using a detailed GEANT-based Monte-Carlo simulation [^] we determine overall 
detection efficiencies of 15 — 46%, as listed in Table |2[ Efficiencies contain secondary 
branching fractions for K° — > Kg — ► TT + ir~ and 7r° — ► 77 as well as rf and u decay modes 
where applicable. We estimate a systematic error on the efficiency using independent 
data samples. 

In Table || we summarize, for each mode, cuts, efficiencies, and the total number of 
events which pass the cuts and enter the likelihood fit described in the next paragraph. 

To extract signal and background yields we perform unbinned maximum-likelihood 
(ML) fits using AE, M, J 7 , \ cos 6b\ (if not used in J 7 ), and dE/dx (where applicable), 
daughter resonance mass (where applicable), and helicity angle in the daughter decay 
(where applicable). The free parameters to be fitted are the asymmetry ((/ — /)/(/ + 
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TV 
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PI 




M B 


5.2-5.3 


5.2-5.3 


5.2-5.3 


5.2-5.3 


5.2 5.3 


5.2-5.3 


|A£| 


< 0.3 


< 0.2 


< 0.2 


< 0.2 


< 0.2 


< 0.2 


|cos6> sph 


< 0.8 


< 0.8 


< 0.9 


< 0.9 


< 0.9 


< 0.8 


Efficiency 


0.40 


0.15 


0.26 


0.27 


0.29 


0.46 


Events in Fit 


6991 


1558 


21337 


395 


10284 


5407 



Table 2: Summary of cuts and number of events used in each mode. 



/)) and the sum (/ + /) in both signal and background. In most cases there is more 
than one possible signal component and its corresponding background component, as for 
instance we fit simultaneously for K + tt° and 7r + 7r° to ensure proper handling of the Kir 
identification information. The probability distribution functions (PDFs) describing 
the distribution of events in each variable are parametrized by simple forms (gaussians, 
polynomials, etc.) whose parameter values are determined in separate studies. For 
signal PDF shapes the parameter determinations are made by fitting signal Monte Carlo 
events. Backgrounds in these analyses are dominated by continuum ee — ► qq events, and 
we determine parameters of the background PDFs by fitting data taken below the T(45) 
resonance or data taken on resonance but lying in the sidebands of the signal region. 
The uncertainties associated with such fits are used later to assess the final systematic 
error. 

5 Results 

(Kir ) In the mode B ± -> i^vr we find a total of 45.6 jlio | even ^ s with an asymmetry 
of ^cp(^vt°) = -0.27 ±0.23. This corresponds to 28.9 ±7.5 K + ir° and 16.8 ±7.5 K~7r 
events. (Here and elsewhere we will quote the yields S and S for the convenience of 
the reader. These are derivative quantities as the fit directly extracts Aqp and S + S.) 
We note that -kit which is analyzed simultaneously but does not have a sufficiently 
significant yield to measure a branching fraction shows an asymmetry of Acp(^"^°) = 
—0.03 ± 0.40. Fig. |l| shows the likelihood function dependence on Acp(Kit ). As a 
cross check we measure the asymmetry of the background events, finding 0.023 ± 0.026 
for Kit background and 0.000 ± 0.017 for 7T7T background. These values are consistent 
with the expected null result for continuum background. 

(K s n) In the mode — ► Kgir^ we find a total of 25.2^'g events with an asymmetry 
of A C p(K s ir) = 0.17 ± 0.24. This corresponds to 10.2 ± 4.0 K s ir + and 14.5 ± 4.4 K s n~ 
events. The background events show an asymmetry of —0.02 ±0.04, consistent with zero 
as expected. 
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30 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




-1.0 -0.5 0.0 0.5 1.0 



Acp 

Figure 1: B — > K ± n° mode: — 2 \n(C/C max ) as a function of *4cp- 



[Kit) In the mode B — > K^tt 1 ^ we find a total of 80.2^^'q events with an asymmetry of 
Acp(Ktt) = -0.04±0.I6. This corresponds to 41.6j|;{j K+n* and 38.6±|;? K~tt + events. 
The dependence of the likelihood function on Acp is shown in Fig. The background 
events show an asymmetry of —0.02 ± 0.04, consistent with zero as expected. 

(co>7r) In the mode — > ujtt^ we find a total of 28. 5^73 events with Acp(utt) = 
— 0.34lg;26. This corresponds to ID.lt^g ujtt + events and 9.4^0 events. Figure ||] 
shows the likelihood as a function of A&P- The fit also allows for a possible charge 
asymmetry in the continuum background, and finds values consistent with zero: — 0.01 3± 
0.015 for ujK ± background, and -0.001 ± 0.01 for uj-k background. 

(i]'K) We find Aqp^'K) = 0.03±0.12. Separating the rj'K signal sample by submodes 
7/ — > r]ir + -K~ and rj — > pj we find.Acp = 0.06±0.17 and .Acp = —0. Of ±0.17, respectively. 
Fig. |I] shows the dependence of the fitted likelihood function on Acp for the separate 
and combined submodes. Background rj'K events are found to have an asymmetry of 
—0.01 ± 0.07 in the r]Tr + ir~ mode, and —0.009 ± O.Of 5 in the pj mode, both consistent 
with zero as expected. 
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-1.0 -0.5 0.0 0.5 1.0 

Acp 

Figure 4: rj'K ± mode: —2 In (C/C max ) plotted versus Acv- The dashed curve shows 
the fit result for the rj' — > rj-K-K decay, the dotted curve for rj' — > p7 and the solid 
curve is the combined result. 

6 Systematic Errors 

The charge asymmetries measured in this analysis hinge primarily on the properties 
of high momentum tracks. The charged meson that tags the parent b/b flavor has 
momentum in all cases between 2.3 and 2.8 GeV/c. In independent studies using very 
large samples of high momentum tracks we have searched for and set stringent limits on 
the extent of possible charge-correlated bias in the CLEO detector and analysis chain 
for tracks in the 2 — 3GeV range. Based on a sample of 8 million tracks, we find Acp 
bias introduced by differences in reconstruction efficiencies for positive and negative high 
momentum tracks passing the same track quality requirements as are imposed in this 
analysis is less than ±0.002. For K^n^ combinations where differential charge-correlated 
efficiencies must also be considered in correlation with K/tt flavor, we use 37,000 D° — > 
Kir(ir°) decays and set a corresponding limit on .Acp bias at ±0.005. These D° decays, 
together with an additional 24,000 decays, are also used to set a tight upper limit of 
0.4 MeV/c on any charge-correlated or charge-strangeness-correlated bias in momentum 
measurement. The resulting limit on Acp bias from this source is ±0.002. We conclude 
that there is no significant Acp bias introduced by track reconstruction or selection. 
We note that for each mode we crosscheck the asymmetry of the background events 
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(normally a fairly large sample) and find results consistent with zero as anticipated. 

Particle identification information for and tt^ is shown in Fig. [|. No significant 
differences are seen between different charge species. Quantification of the effect on *4cp 
is covered by PDF variation studies discussed immediately below. 





-5 -2.5 2.5 5 -5 -2.5 2.5 5 

(dEdx — expected)/a (dEdx — expected)/a 



Figure 5: Normal central dE/dx distributions for and tt ± . 



All PDF shapes which are used in the maximum likelihood fits are varied within lim- 
its prescribed by the fits which determine the shape parameters to assess the systematic 
error associated with uncertainty in the parameters. The resulting changes in *4cp are 
summed in quadrature to estimate the systematic error due to possible misparametriza- 
tion of the PDF shapes. This contribution to the systematic error ranges from ±0.02 
to ±0.04 depending on mode. 

We choose to assign a conservative systematic error ±0.05 for all modes. 



7 Summary 

Table ^ and Fig. ^ summarize all results. 

We have presented in this paper the first measurements of charge asymmetries in 
hadronic B decay. We see no evidence for CP violation in the five modes analyzed here 
and set 90% CL intervals (systematics included) that reduce the possible range of .Acp 
by as much as a factor of four. While the sensitivity is not yet sufficient to address the 
rather small -4cp values predicted by factorization models, extremely large Acp values 
which could arise if large phases were available from final state interactions are firmly 
ruled out. For the case of Kir and rj'K we can exclude at 90% confidence values of |-4cp| 
greater than 0.35 and 0.28 respectively. 
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Mode 


s 


S 




90% CL 
interval 


Ktt° 


16.8 ±7.5 


28.9 ±7.5 


-0.27 ±0.23 ±0.05 


[-0.70,0.16] 


K° s ir 


14.5 ±4.4 


10.2 ±4.0 


0.17 ±0.24 ±0.05 


[-0.27,0.61] 


Kn 


38.6i|;S 


41.6i£§ 


-0.04 ±0.16 ±0.05 


[-0.35,0.27] 


rj'K 


51.7 ±9.2 


48.7 ±8.9 


0.03 ±0.12 ±0.05 


[-0.22,0.28] 


UJ7V 


9.4£tg 


i9.i±S5 


-0.34 ±0.25 ±0.05 


[-0.80,0.12] 



Table 3: Acp measurements in five charmless B decay modes. CLEO 1999 Pre- 
liminary. 



The search for CP violating asymmetries will intensify with the new datasets expected 
in the coming years. The precision of such searches will be statistics limited and should 
improve with integrated luminosity as 1/y/C- If modes with large asymmetries and 
reasonable branching ratios exist they could be found within a few years. 
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